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^ ' Abstract In this paper, we give three equivalent properties of the class of multivariate 

+^ ■ simple cubic natural exponential families (NEF's). The first property says that the cu- 

mulant function of any basis of the family is a solution of some Monge- Ampere equation, 
the second is that the variance function satisfies a differential equation, and the third is 
characterized by the equality between two families of prior distributions related to the 
J> ' NEF. These properties represent the extensions to this class of the properties stated in [1] 

and satisfied by the Wishart and the simple quadratic NEF's. We also show that in the 
. real case, each of these properties provides a new characterization of the Letac-Mora class 

; of real cubic NEF's. 

od ■ 

■ Keywords: Natural exponential family, variance function, cumulant function, Monge- 
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5^ : 1 Introduction and preliminaries 

For the convenance of the reader, we first introduce some notations and recall some facts 
concerning the natural exponential families and their variance functions, our notations are 
the ones used in [10]. Let ii^ be a linear vector space with finite dimension n, denote by 
E* its dual, and let E* x E M. : {9,x) i-^> {6,x) be the duality bracket. 
If /i is a positive radon measure on E, then 

L^{e) = [ exp{{e,x))fiidx) (1.1) 

JE 

denotes its Laplace transform. We also denote by ^A{E) the set of measures fj, such that 
the set 

e{n) = interiorie* € E*; L^{e) < +00} (1.2) 
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is non empty and fi is not concentrated on an affine hyperplane of E. The cumulant 
function of an element /j, of A4{E) is the function defined for 6 in by 

ki,{9) = log L^{e). 

To each /j, in A4{E) and 9 in Q{p), we associate the probability distribution on E 

P{e,lj){dx) = exp{{e,x) - k^{9)) fiidx). 

The set 

F = F{fi) = {p{e,fiy, eeeifi)} 

is called the natural exponential family (NEF) generated by fi. We also say that is a 
basis of F. 

The function is strictly convex and real analytic. Its first derivative fc^ defines a 

diffeomorphism between 6(/x) and its image Mp- Since A:^(^) = J xP{6, n){dx), Mp is 

called the domain of the means of F. The inverse function of k'^ is denoted by ■0^ and 
setting 

P{m,F) = P{^l;^{m),fi) 
the probability of F with mean m, we have 

F = {P{m,Fy, m G Mp} , 

which is the parametrization of F by the mean. 

Now the covariance operator of P{m,F) is denoted by Vpim) and the map 

Mf Ls{E\E)- m ^ VF{m) = k'l^i^m)) 
is called the variance function of the NEF F. It is easy proved that for all m G Mp, 

VF{m) = {i^'^{m))-\ 

and an important feature of Vp is that it characterizes F in the following sense: If F 
and F' are two NEFs such that Vpim) and VF'{m) coincide on a nonempty open set of 
Mp n Mp', then F = F' . 

Now, let us examine the influence of an affine transformation and a power convolution 
on a NEF F = F{ji). If ip{x) = a{x) + b, where a G GL{E) and h e E, \s an affine 
transformation of E, then ip{F{^)) = F{ip{fi)), M^(^p^ = (p{Mf), and 

V^(F)i'm) = a VF{ip~^{m)) a*, 

where a* is the transpose of a. On the other hand the set 

A(/x) = {A > 0; 3 /XA e M{E) such that L^^{e) = {L^{e)f for ah 9 G e(/x)} 

is called the Jorgensen set of jx and the measure ^\ is the A— power of convolution of /x. 
For A in A(/x), we have that 

Mf{^^) = XMf, and VfJ™) = A VFi^)- 



2 



A very interesting fact is that the most common real and multivariate probability 
distributions belong to the natural exponential families such that the variance function 
is a polynomial of degree less then or equal to three in the mean m. For instance, up to 
affine transformations and power of convolution (up to the type), the Gaussian, Poisson, 
gamma, binomial, negative binomial and hyperbolic cosine distributions form the class 
of all real NEF's whose variance function is a polynomial of degree less than or equal 
to 2 characterized by Morris [14]. Letac and Mora [11] have added six others types of 
distributions, namely, the inverse Gaussian, Ressel, Abel, Tackas, strict arcsine and large 
arcsine, to get the class of real cubic NEF's, that is the class of NEF's such that variance 
function is a polynomial of degree less than or equal to three. The classification of NEF's 
with polynomial variance function have been extended to the multivariate NEF's. The 
multivariate version of the Morris class, called the class of simple quadratic NEF's, has 
been completely described by Casalis [1], it contains 2n + 4 types. Hassairi [6] has defined 
and characterized the so-called class of multivariate simple cubic NEF's which is the 
natural extension of the class of real cubic NEF's. It is worth mentioning here that the 
simple quadratic NEF's are not the only families which have quadratic variance functions, 
the Wishart families on symmetric matrices have also quadratic variance functions. The 
classifications of NEF's by the form of the variance function provide an important tool in 
the study of distributions. In fact, in many important cases, the variance function is very 
simple and is easier to use than the distribution itself or the Laplace transform. Moreover, 
the fact that the variance function is quadratic or cubic, is not only a question of form, 
but the form corresponds to some very interesting analytical characteristic properties. In 
this respect, let us mention that for the Morris class of real quadratic NEF's, we have the 
Meixner characterization based on some families of orthogonal polynomials which generate 
exactly the Morris class (scc[12]). Another characterization due to Fdnsilver[4] states that 
a certain class of polynomials naturally associated to a NEF is orthogonal if and only if the 
family is in the Morris class. This characterization has been extended to the Casalis class of 
simple quadratic NEF's by Labeye-Voisin, and Pommeret [9]. Concerning the cubic NEF's, 
Hassairi and Zarai [7] introduced a notion of 2-orthogonality for a sequence of polynomials 
to give an extended version of the Meixner and Feinsilver characterization which subsume 
the Letac-Mora class of real cubic NEF's. Hassairi and Zarai [8] have also introduced a 
notion of trans-orthogonality for a sequence of multivariate polynomials to extend their 
characterization result to the class of multivariate simple cubic NEF's. Besides these 
characterizations based essentially on different notions of orthogonality of polynomials, 
it is stated in Casalis[l] that the simple quadratic NEF satisfies a property based on 
two conjugates families of prior distributions related to the NEF. For a NEF F = F{^), 
consider the family of prior distributions H introduced by Diaconis and Ylvisaker [3] and 
defined by 

n = {7rt,^o(d^) = Ct,mo expt((mo,^) - k^{e)) le^^){e)de,t > 0,mo G Mp} (1.3) 

where Ct,mo is a normalizing constant. Consider also the family H* introduced by Consonni 
et al [2], see also [5] and defined by 

n* = K™o>*>o>"^oeM^} (1.4) 

where 

T^t,mo{dm) = Cl^^ expi( {mQ,tpi^{m)) - kij,{tpi^{m)) ) lMp{m)dm, 
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and the constant C*,„jp is a normalizing constant. Then, when F{p,) is a Wishart or a 
simple quadratic NEF, we have that A:^(n) = 11*. It is also shown that this property is 
equivalent to two other properties expressed in terms of some differential equations satisfied 
by the cumulant function k^. In the real case, the property characterizes the Morris class 
of real quadratic NEF's, that is k'^{Il) = 11* if and only if the NEF F is in the Morris 
class. In the present paper, we extend these results to the class of multivariate simple 
cubic NEF's. We construct two families of prior distributions related to a multivariate 
NEF, and we show that these families coincide when the NEF is simple cubic. We then 
show that this property is equivalent to the fact that the cumulant function is a solution 
of some Monge- Ampere equation and also equivalent to the fact that the variance function 
satisfies a differential equation. As a corollary, we obtain three new characterizations of 
the Letac-Mora class of real cubic NEF's. 

2 Some equivalent properties 

Throughout this section, we suppose that F = F{p,) is a NEF on a linear vector space E 
with dimension n. Besides the family 11 of prior distributions defined in (1.3), we introduce 
another family n of prior distributions. Let P be in E* such that the set 

G = {eeeif,y,i + {p,k'^{e))>o} 

is nonempty, and denote M = fe^(S). Consider the family of prior distributions 

n = { nt,mo ; t e m; , mo G Mp}, (2.5) 

where 

T^t,mo{dm) = Ct,mo (1 + (^,m))~"~^ ex.pt{{mo,tpfi{m)) - A;^('0;u(m))} l^{m)dm. 
With these notations, we next state and prove our first main result. 

Theorem 2.1 The three following properties are equivalent 

(1) There exists (a, b,c) E: E x such that for all m in Mp, 

detyF(m) = (1 + {P,m)y+^ e-x^{{a,^p^,{m)) + bk^{^p^,{m)) + c). 

(2) There exists {a,b) E E x M such that for all m in Mp and any basis, (cj)"^^ of E, 
with dual basis (e|)"^i, we have 

j:mm){eiM = i + + ^ + bm. (2.6) 

(3) There exists an open subset ft of M*^ x Mp such that 

fe^(n) = flo = {Tft,™^; (i,mo) G O}. 
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Note that (1) may be stated in terms of the cumulant function as there exists (a, 6, c) G 
E xJR^ such that for all 9 in e(//), 

detA:;:(^)= {l + {P,k'^ie))r+^ eM{a,0) + bk^ie) + c}, 

that is the cumulant function is solution of some Monge- Ampere equation (see [15]). 
Proof We will show that (1) <^ (2) and (1) <^ (3). 

(1) =^ (2) Suppose that Vpim) satisfies (1), then we have 

log det VfM = (n + 2) log(l + m)) + {{tpi,{m),a) + b A;^(V'^(m)) + c}. 

Taking the derivative, we get 

which is equivalent to 

jZ[yF{m){.)V^\m){e,)]{el) = + + {^l^'^{m){.),a) +h{m,^l^'^{m){.)). 

Replacing (.) by VF'(m)(.), and using the condition of symmetry 

V\m){V{m){a)){P) =V'{m){V{m){P)){a) Va,/5 G E* (see [10],pa5e 103), 
we obtain 

t[v;(»)(e.)(.)l(e-) = + + (a. (.)) + Km, (.))• (2.7) 

As Vpim) is symmetric, we get 

E[^;M(e,)](e:) = .^'!.^\ v>(m)(/3)+a + bm. 

(2) ^ (1) Suppose that (2) holds. Then, we easily get (2.7). 
Replacing, in (2.7), (.) by Vp^{m){.), one obtains 

fy^{m){ei)Vf\rn){.){el) = + (a, ^;(m)(.)) + 6(m, ^;(m)(.)). 

This is equivalent to 

trace{Vf\m)V'{m){.)) = ^Y+l/^fm)"^^ ^ {a,^',{m){.)) + b{m,^|;'^{m){.)). 
Integrating, we deduce that there exists c in IR such that 

logdet(VF(m)) = (n + 2) log(l + (^, m)) + {(V'M(m), a) + b k^{il^t,{m)) + c}, 
and the result follows. 
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(1) =4> (3) Suppose that (1) holds, and define 

n = {{t, mo) G ]R+ X Mf ; t > 6 and uiq € {I - ^)Mf - 

Take {t,mo) in J7 and denote u the image of ^t,mo by V'/u- Then it is easy to verify that 

i^ide) = Ct,m^ exp{(imo + a,^) - (t - 6)A;^(0)} \^{&)dB. 

Since (t, mo) is in wc have that t — 6 > and ^ g Mp. 

t — 

tmo + OS 1 . 1 

Thus taking = t — o and m\ = ; — , we obtain that 

t — o 

v{dQ) = Ct,,m, exp{t((mi,^) - k^{e))} l^{9)d0. 
Hence 4^fj,{Iln) C IT, and it follows that 11^ C k'^{Il). 

Conversely, if irt^mo is an element of H, then its image a by k'^ is given by 

a{dm) = Ct,mo e"'' {1 + {P^m))"^-'^ eyi^{{tmQ- a,il)^{m)) - {t + h)k^i(il)^{rn))} l~(m)dm. 

Taking t\ = t + h and mi = ^ . Then (ii, mi) is in f2, and we have 

a{dm) = Ct^^mi (1 + (Z?,?"))"""^ expti{(mi,V'^(m)) - fe^(V'^(m))}lj^(m)(im , 

which is an element of Uq . 

(3) => (1) Suppose that k'^{H) = Hq. Then, for an element TTt^mo of H, we have on the 
one hand, 

k'i^{'^t,mo)idm) = (det VF(m))"^ Ct,mo exp i{(mo, V"//!™)) - ^^(V'/iC™))} 'i-^{m)dm. 
On the other hand, since k'^{nt,mo) is in Hq, there exists (ti,mi) in U such that 

k'i^in,mo){dm) = Ctj_,mi (1 + (/3,to))"""^ expti{(mi,V'/,(TO)) - fc^ (m))} lj;j(m)dm. 
Comparing these two expressions of A;^(7rt,^o) gives 

dety(m) = (l + (/3,m))"+2 exp{(a, ^^(m)) + 6fc^(V'^(m)) + c}, 
where a = tmo — timi, b = t\ —t, and c = log( _*''"° ). 

3 Characterizations of the Letac-Mora class of 
real cubic NEFs 

In this section, we prove that a multivariate simple cubic NEF satisfies the properties in 
Theorem (2.1), and that the real versions of these properties characterize the real cubic 



6 



NEFs. Recall that a simple cubic NEF is obtained form a simple quadratic NEF by the 
so-called action of the linear group GL(]R x E) on the NEFs of E. For more details, we 
refer the reader to [6], where a complete description of this class is given. This action 
is in fact an extension of the way in which the Letac-Mora class of real cubic NEFs is 
obtained from the Morris class of real quadratic NEF's. For our purposes here, we need 
only to mention that, up to affine transformations and power of convolution, a simple 
cubic variance function is of the form 



V{m) = {l + {P,m)) (I + m^P) Vi{- 



m 



(3.8) 



'l + (^,m)' 

where Vi is the variance function of a simple quadratic NEF Fi, and m is in {Mp-^^)fs, where 



m 



{Mf,)i3 = {m e Mf,;1 + {P,m) > and - 

l + {p,m} 

The relation (3.8) is invertible and conversely, we have 



Vi(M) = (1 - (^, M)) {I-Mi^P)V{ 



M 



1-{P,M) 



){I-P®M), 



(3.9) 



where M is in (Mi?)_^. 

We also mention that the relation between a simple cubic NEF F{^) and a simple quadratic 
NEF F{v) may also be expressed in terms of the cumulant functions by 



or equivalently by 



= K{e) - ko 

A = -PK{9) + e-Xo 

K{e) = k^{\) - ki 



(3.10) 



(3.11) 

_ 9 = pk^{\) + \-9^ 

where (fcojAo) and (/ci,^i) are constants in M x Note that if /3 = in (3.8) wc obtain 
the simple quadratic class. Then for more accuracy we exclude this case and we keep only 
^mE*\ {0}. 

We now prove that the multivariate simple cubic NEF's satisfy the properties in Theorem 
2.1. 

Proposition 3.1 Let F = F{ijl) he a simple cubic NEF on E, then there exists (a, 6, c) 
in E* X such that 

det(yp(m)) = (1 + (/3,m))"+2 exp{(V';.(m), a) + b k^{i^^{m)) + c}. 

Proof Given that the family F is simple cubic, then there exist /3 in E* and Fi = F(y) a 
simple quadratic NEF such that 



VF{m) = (1 + m)) {I + m®^) Vf,{ 



m 



l + (/3,m) 



{I + I3®m), 
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see (3.8). As det(/ + m /3) = 1 + (/?, m), we obtain 

det(V>(m)) = (1 + (/3,m))"+2 det(y^, ( ^ ^ )). 

Wc now use the fact for a simple quadratic NEF Fi (see [1]), there exist a' in and 6', c' 
in ]R such that, for all M in Mi;'^, 

detiVpAM)) = exp{(a', V.(M)) + h'K{^,{M)) + c'}. 

It follows that 

det(Ff (m)) = (1 + (;3,m))"+2 exp{(V',(^-p^^), a') + h' K{My^^^)) + ^'i- 

777- 

Prom (3.10), putting A = '^^{m) and 9 = ipi,{ j- r- ), we get 

1 "I" \Pi '^1 

Tfl 

^ ^ )) = kf,{tp^{m)) + ko, 

TTl 

i^v{ ^■^ ) = ^'^("i) + /3fe;,(V'M("i)) + Pko + Ao. 

Then 

detyF(m) = (l + (Am))"+2exp{(V';.(m),a') + (6' + (a',/?))/c^(V^(m)) 
+{b'ko + {a',l3ko + \o) + c')}. 

Setting a = a', b = b' + (a',/3) and c = b'ko + {a',/3ko + Aq) + c', we obtain the desired 
result. □ 

As the Letac-Mora class of real cubic NEFs is nothing but the simple cTibic class, 
when the dimension n is equal to 1, this class satisfies the real version of the properties in 
Theorem 2.1. We will show that, in this case, these properties are characteristic. 

Theorem 3.2 Let F = F{fj,) be a NEF on the real line, then F is cubic if and only if 

k'^in) = fx. 

Proof Suppose that fc^(n) = n. Then according to Theorem(2.1), the variance function 
Vpim) satisfies the differential equation 

(1 + l3m)Vp{m) - 3/3 Vpim) = (a + bm){l + /3m). 

Solving this differential equation by standard methods gives 

Vpim) = A (1 + /3m)3 - ^ (1 + Pmf + (1 + /3m), 

which is a polynomial of degree less then or equal to 3. □ 
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